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ABSTRACT 

The old suggestive observation that black holes often resemble lumps of fluid has recently been 
taken beyond the level of an analogy to a precise duality. We investigate aspects of this duality, and 
in particular clarify the relation between area minimization of the fluid vs. area maximization of 
the black hole horizon, and the connection between surface tension and curvature of the fluid, and 
surface gravity of the black hole. We also argue that the Rayleigh-Plateau instability in a fluid tube 
is the holographic dual of the Gregory-Laflamme instability of a black string. Associated with this 
fluid instability there is a rich variety of phases of fluid solutions that we study in detail, including 
in particular the effects of rotation. We compare them against the known results for asymptotically 
flat black holes finding remarkable agreement. Furthermore, we use our fluid results to discuss the 
unknown features of the gravitational system. Finally, we make some observations that suggest 
that asymptotically flat black holes may admit a fluid description in the limit of large number of 
dimensions. 
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1 Introduction 



The tantalizing similarities between, on the one hand, the physics of black holes, and on the other 
hand, the properties of soap bubbles and fluid lumps, have been observed from old and indeed 
motivated approaches like the membrane paradigm. A black hole is bounded by a smooth horizon 
with uniform surface gravity, in analogy with the surface of a fluid droplet which, under the effect of 
surface tension, assumes a smooth shape with constant mean curvature. The negative specific heat 
of a black hole, and its concomitant Hawking radiation, bear a resemblance to the evaporation of 
a fluid droplet. More recently, refs. [UEllS] pointed out a striking similarity between the Rayleigh- 
Plateau instability of a column of fluid and the Gregory-Laflamme instability of black strings [H [5] . 

However, when one tries to take these observations beyond the level of mere analogies, a number 
of apparent difficulties arise. A soap bubble, or a lump of fluid — at least a static one — strives 
to minimize its surface area, whereas a black hole tends to maximize the area of its horizon. It is 
also difficult to make precise in which sense the surface gravity relates to the surface tension of the 
fluid, or whether it should be related instead to the mean curvature of the fluid surface. Underlying 
these difficulties is the question of whether there is any regime in which the gravitational dynamics 
of black hole horizons can be precisely mapped into the dynamics of fluids. 

The issue can be decisively settled within the context of the AdS/CFT correspondence, which 
turns the analogy into a precise duality at least for a certain class of black holes in Anti-de Sitter 
space. Refs. [6l[7] have shown that large black holes in AdS gravity can be mapped to solutions of 
the Navier-Stokes equations of a conformal fluid. In this framework, large black holes in a space- 
time asymptotic to AdS with boundary B correspond to fluid configurations filling the spacetime 
B. However, the physics of these black holes differs significantly from the physics of, say, the 
Schwarzschild black hole in asymptotically flat space. For instance, large AdS black holes have 
positive specific heat, and thus do not disappear via Hawking emission. In fact, the analogies 
mentioned above seem to lose some of their motivation since the fluid duals to these black holes 
do not have a bounding surface. Also, black strings exist in AdS^, with horizon topology M x S'^~'^ 
extended in a direction along an asymptotic boundary of topology Mt x M x 5"^~^ [H [U [10], but 
when these black strings are large they do not suffer from the Gregory-Laflamme instability [TT] . In 
dual terms, they correspond to a fluid filling the boundary. This fiuid is locally stable under plasma 
oscillations {i.e., sound waves propagating in the fluid) so in the long- wavelength approximation 
these configurations are stable. Thus, these set-ups do not seem adequate to make precise the 
connection between black holes and fiuid lumps, and in particular the instabilities that may beset 
them. We need the dual fluid to admit a bounding surface. 

Such bounded fiuid configurations arise naturally for gauge theories with a deconflnement tran- 
sition of first order. At temperatures just above the deconfinement transition, lumps of deconfined 
plasma, described at long wavelengths as fiuids, can exist separated by a domain wall from the 
confining vacuum [12]. A gravitational dual of a theory with this behavior is provided by Scherk- 
Schwarz (SS) compactified AdSd+2 gravity, whose solutions asymptote to M'^"^ x S^, with 
standing for the distinguished SS circle. Under the approximation of a very thin domain wall, 
there are solutions for finite lumps of deconfined plasma whose dual gravitational solution must 
correspond to black objects localized near the infrared end of SS-AdS. These solutions are perfectly 
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regular, provided that the SS circle shrinks to zero size in the bulk, in correspondence with the 
domain wall. This implies that the topology of the corresponding event horizon is given by the 
fibration of the SS circle over the plasma lump geometry, with the circle shrinking to a point on 
the boundary. Plasma balls and plasma rings correspond, respectively, to spherical black holes and 
to black rings in the bulk. These were studied in [13] with an emphasis on the AdSs case, and the 
full phase diagram for balls and rings in AdSg was obtained more recently in [l4j. Crucially, the 
horizon of the black hole is not mapped to the boundary of the fluid, but to its entire bulk. Then 
the apparent contradiction between minimal fluid surface area and maximal horizon area can be 
resolved. 

Thus, this seems to be the correct framework to address the correspondence between black holes 
and lumps of fluid. One aim of this paper is to clarify the puzzles discussed above pertaining the 
variational principles that determine the equilibrium configurations of a fluid and of a black hole, 
as well as the relation between physical magnitudes in both cases. 

Having dealt with this, we attempt to uncover in the physics of these plasma lumps some of 
the characteristic phenomena that higher-dimensional black holes exhibit (see e.g., \15 \ \16 \ \T7\). In 
particular, we focus on plasma tubes, which are dual to black strings. As mentioned above, fluid 
analogues of the Gregory-Laflamme instability and of non-uniform black strings [18^ [T9] that branch 
off at its threshold have been investigated recently [U O O [20l Ell [22] • Here we reproduce and 
extend these results to include rotation, which uncovers novel phases of non-uniform tubes. But 
now we can regard these results as more than an analogy: following [7] they are actual solutions 
to a controlled approximation to the Einstein equations in SS-AdS and thus describe gravitational 
dynamics of black holes and black strings localized near the infrared in the SS-AdS spacetim^. 

It turns out that these results show remarkable agreement with those found for black holes 
and black strings in asymptotically flat space. This prompts the question of whether the latter 
should admit, in some limit, a precise description in terms of fluid dynamics. After all, even if 
the equivalence between AdS black holes and fluids is very satisfying and promising, it does not 
yet explain the observations that suggested an analogy between asymptotically flat black holes 
and fluids in the first place. While we do not have a complete answer to this, we will make some 
observations that suggest that asymptotically flat black holes may be accurately described by fluid 
equations in the limit of very large number of dimensions. 

The plan of the rest of the paper is the following: Section [2] reviews the basics of the relativistic 
hydrodynamics of plasma lumps, and provides the dual identification between temperature of the 
black hole and parameters at the surface of the fiuid. Section [3] discusses the equivalence between 
several variational principles for fluids. This allows to understand how maximization of the entropy 
of a black hole can be equivalent to minimization of the surface area of a fluid. In section [5] we 
explore static and rotating configurations of plasmas in equilibrium in a space with one (large) 
compact direction, with a particular emphasis on tubes extended along this direction. In sectionlU a 
detailed analysis of the Rayleigh-Plateau instability for these is performed. We obtain the dispersion 
relations for the unstable modes, both for static (Fig. [5]) and spinning tubes (Fig. [6]). After these 

^Strictly speaking, the analysis of [7] refers to black branes dual to fluids that fill the spacetime they live in. We 
can reasonably expect that the extension to the situation of interest here should only be technically more complicated 
due to the absence of analytic solutions for plasma balls. 



3 



analyses, we discuss in section [6] their range of validity, and in section [7] how they compare to 
the phases and stability of black strings in vacuum gravity. Section [8] revisits the issues posed at 
the beginning of the paper and suggests that vacuum black holes in the limit of large number of 
dimensions may admit a fluid description. 

Note added: during the final stage of this project we have been informed of the work [23j, which 
overlaps with some of our results. 



2 Hydrodynamic description of deconfined plasma lumps 
2.1 Navier-Stokes and Young-Laplace equations 

Large black holes in AdS correspond to deconfined gluon plasma in the boundary CFT. In the long- 
wavelength approximation, this plasma behaves as a fluid, whose properties can be deduced from 
its gravitational bulk description. In this approximation, one can perform a derivative expansion, 
whose leading order is given by a perfect fluid, while next orders introduce dissipative and diffusive 
phenomena. More precisely, it has been shown that the long-wavelength sector of pure AdS gravity 
is described by Navier-Stokes equations on the boundary, and the stress tensor of the corresponding 
fluid has been computed up to second order in the derivative expansion in five [21 [2^ and four 
dimensions [25j . 

We are interested in the case where the boundary theory admits a confining phase, which close 
to the deconfinement temperature is separated from the deconfined plasma by a thin domain wall 
|12j . We assume that our plasma lumps are large enough to neglect the thickness of the wall. The 
analysis of this subsection follows largely [13] . We shall neglect subleading dissipation and diffusion 
contribution^. In this approximation, the deconfined plasma fluid behaves as a perfect fluid, and 
the leading order stress tensor is given by the sum of a perfect fluid part and a boundary surface 
contribution, describing the capillarity of the fluid, 

^ ~ ^perf "I" -'bdry ' 

T;:.i = [{P + P) ^'^^ + Pan ©(-/) > TZry = -<^hndf\ 6{f) . (2.1) 

Here, is the fluid velocity, p, P and a are its density, pressure and surface tension respectively. 
The fluid boundary is defined by f{x^) = 0, has unit normal = d^f /\df \ pointing towards the 
confining phase, and h^'^ = g^" — n^n'^ is the projector onto the boundary. The velocity field is 
subject to the boundary condition 

w''n^ = 0, (2.2) 

requiring that the fluid velocity must be orthogonal to the boundary normal or else the fluid would 
not be confined inside the boundary but flow through it. Finally, G(— /) is the Heaviside function 
(0(— /) = 1 inside the fluid and zero everywhere else). 



^In particular, for stationary equilibrium solutions, these contributions vanish. We will comment on the effects of 
higher order corrections on the analysis of the perturbations near the end of the article. 
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The equations describing the dynamics of the fluid fohow from the conservation of the energy 
momentum tensor p.ip . The volume and boundary contributions must vanish independently; the 
former, projected along the fluid velocity Ui, yields the relativistic continuity equation, 

7x^V^/5 + (p + P)V^u^ = , (2.3) 

(where V ^u^ = 6 is the expansion of the fluid) that we use to reduce the remaining equations to 
the relativistic Navier-Stokes equation, 

(p + P)u''VX = - ia^" + ^^^") • (2.4) 

Here u^V ^u'^ = is the acceleration of the fluid and g^^ + u^^vy is the projector onto the subspace 
orthogonal to u^^. 

Now consider the boundary contribution to the conservation of the stress tensor, 

[{p + P)u^u'' + \dfK + '^V^ (h^dfl) = . (2.5) 

Projecting this equation onto n^, using (j2.2p . and integrating across the boundary we get 

P<- P>=aK , K = h^^'V^nf' , (2.6) 

where K is the trace of the boundary's extrinsic curvature (twice the mean curvature H of the 
bounding surface), and P< — P> is the pressure jump when we cross the boundary from the exterior 
into the interior. Equation (j2.6p is the Young-Laplace equation that describes the capillary pressure 
difference AP sustained across the interface boundary due to surface tension a |26j . 

2.2 Stationary plasma configurations 

The Young-Laplace equation we found in the last subsection expresses the balance of forces on 
the plasma boundary and holds in the most general dynamical situation. We will turn now to the 
study of stationary plasma configurations in hydrodynamical and thermodynamical equilibrium. 

To begin with, we need to characterize stationary configurations of plasma. Let us assume 
that the background spacetime is stationary, with stationarity timelike Killing vector ^, and a set 
of commuting, linearly independent spacelike Killing vectors xiy corresponding to the isometrics 
of the background. Of the latter, the subset that commutes with the velocity field will also be 
symmetries of the fluid, C^jU = 0. The associated conserved charges will be the energy E and the 
linear /angular momenta Jj of the plasma associated to fluid symmetries. 

Let V be the spatial region filled by the plasma, dV its boundary, and the velocity field of 
the fluid. Stationarity, for an isolated system, requires that there be no dissipation, and therefore 
the bulk and shear viscosity terms of the stress tensor must vanish. This condition is met if and 
only if the expansion 9 and the shear a^'^ of the fluid both vanish. Under these assumptions, one 
can show that there exists a function a such that au^ is a Killing vector. 

Indeed, for a shearless, expansion-free velocity field, 

V^'u^' = uj^'' - u^a\ (2.7) 
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where oj^'^ is the vorticity of the fluid and its acceleration. Moreover, a fluid with local entropy 
density s and local temperature T has to satisfy the Euler relation, 

p + P = Ts. (2.8) 

Differentiating this equation and using the first law of thermodynamics, we get the Gibbs-Duhem 
relation dP = sdT . Hence, since for a stationary configuration u^'SJ = 0, the Navier-Stokes 
equations can be rewritten as 

a^ = -(/9 + P)-iV^P = -V^lnr, (2.9) 

where, in the last step, we used both the Euler and the Gibbs-Duhem relations. Note that this 
equation expresses simply the fact that the heat flux vanishes 

= -K{g^'^ + u^n'^)(V^T + a^T) = (2.10) 

or, in other words, stationary plasma configurations are both at hydro dynamical and thermal 
equilibrium. Substituting equation (j2.9|) in (|2.7p it follows that 

V(^(an^)) = an(^V^) In(aT). (2.11) 

Hence, if we choose the function a to be inversely proportional to the local temperature of the 
plasma, a = T/T with T an integration constant, the vector field au^ solves the Killing equations, 
and must therefore be a linear combination of the Killing vectors ^ and xJfl- Hence, stationary 
configurations have a velocity field given by 

n = ^ (e - ^ixi) ■ (2.12) 

The Killing vectors appearing here are a subset, although not necessarily the full set (since there 
may be rotational symmetries with no rotation velocity of the fiuid in their direction), of the 
abelian symmetries of the fluid, which themselves are a subset of the symmetries of the background. 
The most general solution is determined by the constant parameters T and 0/, since the local 
temperature T is then fixed by the relation = — 1 to be 

T = 7r, (2.13) 

where 

7"' = V-i^-^iXi? (2.14) 

is the redshift factor relating measurements done in the laboratory and comoving frames, and the 
constant T is the equilibrium plasma temperature. It follows from the general form (j2.12p of 
that a stationary equilibrium fluid configuration has to be in rigid roto-translational motion. 

Combining the Euler relation (j2.8p and the Young-Laplace equation (j2.6p . we can relate the 
temperature parameter T of the fiuid (which, for a non-boosted lump in an ultra-static background 

^Actually, it could depend linearly on Killing vectors not commuting with xi- However, one can always choose 
the set of Killing vectors xi to be adapted to the motion of the fluid under consideration. 
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spacetime with = —1, is the temperature of the fluid at the axis of rotation) to a combination 
of several magnitudes at the fluid surface, 

r="^ + ^ (2.15) 

In the duality to a black hole, T corresponds to the Hawking temperature of the horizon. We see 
that T is not simply proportional to the surface tension or to the mean curvature, although it grows 
linearly with either of them. Note also that for a static fluid K will be constant over the surface, 
but in a stationary configuration the curvature K will adjust itself to variations of the fluid velocity 
near the boundary. Let us also mention that the angular velocities of the black hole horizon are 
identified with the 0/. 

So far we have not specified the equation of state of the fluid. We are interested in the d- 
dimensional (non-conformal) plasma describing the hydrodynamic limit of the {d + l)-dimensional 
CFT dual to Scherk-Schwarz reduced AdSrf+2 gravity, whose equation of state in d = n + 2> dimen- 
sions reads (see [13]) 

n + 6 

and 

,.,„,.K/,..«.(^)-. r.(^)*, p..) 

with po and a constants. This equation of state is normalized such that the vacuum pressure 
outside the plasma vanishes. When n = 0, if we consider the theory that is dual to type II strings 
in SS-AdSs x 5^, we have 

The surface tension of the domain wall has been computed numerically in |12j for n = and n = 1. 
In units of Pq/T^ they obtain 

c7„=o = 2.0^, and ^^=1 = 1.7^ . (2.19) 

J-c 

For a plasma in equilibrium, it follows from (j2.13p and (|2.16p that the density and pressure 
satisfy the relations 

p = p.7"+' + Po , P = ^ 7"+' - Po , (2.20) 

where 

= (n + 3)ar"+^ (2.21) 

is the difference p — po between the plasma energy density and the energy density of the confining 
vacuum at points where 7 = 1. We see that all the dependence of the density on the position is 
entirely given by the 7 factors. 
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3 Variational principles for equilibrium plasma configurations 



In this section we will show that the Young-Laplace equation (j2.6p can also be obtained using 
two equivalent variational principles, namely through a maximization of entropy and/or through a 
minimization of potential energy. As an easy byproduct we will see that these are also equivalent 
to minimization of the surface area for static fluids. Our formulation will be relativistic, fully 
covariant, and valid for a large class of stationary background spacetimes. While the relation 
between these two variational principles is probably well-known, we were not able to find it in the 
literature and we surmise that the relativistic covariant generalization that we present is original. 

3.1 Relativistic soap bubbles and second law of thermodynamics 

Equilibrium plasma configurations must have constant temperature and are expected to be those 
that extremize the entropy of the plasma while keeping its energy and momenta fixed. We will use 
this variational principle to give an alternative derivation the Young-Laplace equation, shedding a 
new light on the relation between black holes and their membrane analogues. 

As explained in the previous section, equilibrium configurations on a stationary background 
spacetime are in rigid motion, and have a velocity field of the form ()2.12p . determined by the 
temperature T of the plasma and the boost parameters Again ^'^ is the stationarity Killing 
vector and we choose a time coordinate t such that ^ = 5j. This allows us to foliate the spacetime 
into constant t hypersufaces S^, that we use to define the conserved charges. We call the 
unit normal vector to these hyper surf aces, and we make the additional hypothesis that the Killing 
vectors xi generate isometrics of St, i.e., k ■ xi = 0- Then, given any Killing vector ip^, one can 
define the associated conserved charges 

Q[V] = [ dv T^.kf^r, (3.1) 

where dv is the induced volume measure on S^. 

To perform the maximization of the entropy at constant energy and momenta, we define the 
action 

i[r] = s[r] - mn + cjiJi[v] , (3.2) 

where (3 and cD/ are the Lagrange multipliers associated to our constraints. The total entropy of 
the fluid is the conserved charge associated to the entropy density current su'^, 

S[V] = - [ {k-u)sdv, (3.3) 
Jv 

while the energy and momenta are given by 

E[v]= [ [{p+p){i-u){k-u) + {k-i)p]dv-<j I k^eh^,\dmf)dv, 

JV JT.t 

Ji[P] = [ iP + P)ik- u)ixi ■u)dv-a [ h^uk^x'i\dfm)dv . (3.4) 

JV JT.t 

In addition to the bulk terms, both the energy E and the momenta J/ have a surface contribution 
proportional to the surface tension. The boundary term in E\P] corresponds to the surface tension 
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potential energy, and reduces to a times the area of the boundary of the fluid lump for a configu- 
ration respecting the spacetime symmetries in an ultra-static background. The boundary term in 
Jl\P\ is present only if • ''^ 7^ so the fluid boundary is not invariant under the action of x/j for 
instance when the surface is not axisymmetric. 
Taking the Lagrange multipliers to be given by 

/? = l/r, lD, = /?J7,, (3.5) 

making use of the form (j2.12p for u and eliminating the entropy density s by using the Euler 
relation, the action simplifies to 

m = -ii I {k-OPdv + Pa [ {k-0\dfm)dv. (3.6) 

To extremize this action, we need to find its variation under a deformation of the region V 
occupied by the fluid. Notice that, since by hypothesis the configuration is stationary, the action is 
time independent, and its extrema are not changed if we integrate it over the time t. Then, noting 
that the spacetime volume element is 

dV = -{i-k)dvdt, (3.7) 
it follows that the quantity we have to extremize assumes the form 

[l[V]dt = p[ PQ{-f)dV -f5a f \dmf)dV, (3.8) 

J JM JM 

where the integrals extend to the full spacetime M, and the resulting functional is manifestly 
covariant. Also, note that the last integral is the area of the region of spacetime (world-volume) 
spanned by the boundary of the fluid. To extremize (jS.Sp under deformations of "P, we compute its 
variation under a change 5f of the function / and integrate by parts the last term, 

5f [ I[V] dt = P [ {P- aK) 6f 6{f) dV . (3.9) 
J JM 

We used here the well-known fact that the variation of the area of an hypersurface is precisely given 
by twice the mean curvature, or its extrinsic curvature K. Therefore, requiring that ()3.9p vanishes 
for any deformation 6f of the boundary, we obtain that the variational principle boils down to 

aK = P, (3.10) 

where P is the pressure jump at the boundary of the fluid. This, for non-vanishing external 
pressure, is again the Young-Laplace equation (j2.6p that we derived in the last section using the 
conservation of the energy-momentum tensor. This time we found that it also follows from requiring 
extremization of plasma entropy for fixed energy and momenta. Note that the derivation assumes 
only the stationarity of the configuration and of the background geometry, and is independent of 
the fluid equation of state. Besides, the derivation does not assume any condition on the shape of 
V, in particular it covers non- axisymmetric cases where the geometry of the fluid configuration does 
not share the symmetries of the fluid motion, as in the case of the two-lobed figures of equilibrium. 
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3.2 Equivalent variational principles for spinning lumps 

A universal behavior of fluids is that they always pick boundary configurations that reduce their 
potential energy for a fixed volume. For static solutions, this implies that the area of the fluid 
surface is minimized. For stationary solutions, the potential energy not only has a surface tension 
term but also a centrifugal contribution. We want to verify the equivalence between this potential 
energy minimization (hence a mechanical formulation) problem and the entropy maximization 
problem (a thermodynamic formulation) studied in the previous subsection. 

To find an energetic formulation of the variational principle, we first observe that, using p.l2p . 
the total action (j3.2p assumes the form ()3.6p in which the second term in the right hand side has 
the interpretation of surface potential energy, while the pressure integral has the interpretation of 
centrifugal potential energy, since the radial variation of pressure measures the centrifugal force 
per unit volume of fluid. 

Therefore, we define a new action 

T[V] = U^[V] + U,f[V] - vV[V] , (3.11) 
where U^lV] = crA[dV] is the surface tension potential energy associated to the area 

A[dV] = - [ {k-0\dmf)dv, (3.12) 

the centrifugal potential energy, obtained by integrating the pressure over the plasma volume is 

UcfiV] = [ {k-0{P-Py)dv, (3.13) 
Jv 

with P> an integration constant, and finally the volume V[V] reads 

V[r] =- [ {k-Odv, (3.14) 

and is kept fixed as we perform the variation through the Lagrange multiplier rj. Note that, in 
an ultra-static background, ^[97^] and V[P] reduce to the usual definitions of boundary area and 
volume of the plasma, respectively. Again, one can then easily verify that the Euler-Lagrange equa- 
tion associated with (j3.1ip yields the Young-Laplace equation (j2.6p . and the Lagrange multiplier 
takes the value rj = P>, where -P> is the exterior vacuum pressure. 

The two variational approaches, maximization of entropy and minimization of area, are equiva- 
lent because the actions (13. 2p and (13. lip are proportional to each other. Indeed, starting with the 
action (|3.2p . successive use of ()3.4p . (13. Sp and (j2.8p yields the following relation between actions 
(lOl and dam), 

I[V] = -pT[V] . (3.15) 

Since the two actions are the same up to a negative constant, minimizing the potential energy 
for fixed volume, SI[V] = 0, is naturally equivalent to maximizing the entropy for fixed conserved 
charges, 6I[V] = 0. 
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Finally note that in static configurations, $7/ = 0, the pressure is constant inside the fluid and 
we can write 

I[V] = -I3I[V] = -Pa {A[crP] - KV[V]) , (3.16) 

and when this action is extremized the Lagrange multiplier K is required to be twice the mean 
curvature. Therefore, non-rotating equilibrium fluids have boundary shapes that minimize the 
surface area for fixed volume and their boundaries have constant mean curvature, a well known 
property of fluids. 

Black holes satisfy the variational principle that their entropy is extremized for fixed energy 
and angular momenta — this is essentially the first law. In the duality between SS-AdS black holes 
and fluid lumps, the entropy, energy and spins are identified on both sides, while the temperature 
is mapped according to eqs. (|2.15p and 0/ are the angular velocities. Then our analysis shows that 
maximization of the black hole horizon area is equivalent, for static configurations, to minimization 
of the fluid surface area. When rotation is present, the connection is less simple geometrically, but 
still easily expressed as an extremization of a functional. 

Thus, in this context the proposal in [IJ to model black strings with non-gravitating cylinders of 
fluids can be made exact. There are, however, some differences between the analysis in refs. [H I20j 
and ours: theirs was a study of non-relativistic fluids, whereas our fluids are relativistic. For static 
fluids the differences are smal0, but for non-static or out of equilibrium situations the relativistic 
corrections become crucial. But the more important differences are that, first, we have a definite 
equation of state for the fluid and a specification of the value of the surface tension. This allows 
us to compute thermodynamic quantities of the fiuid that are identified with those for the black 
hole, something that was impossible in [UEO]. Second, the number of dimensions for the fiuid and 
the black hole are different: a fluid in d spacetime dimensions is mapped onto a black hole in 2 
dimensions. In particular, the entire volume of the fluid can be regarded as mapped onto the black 
hole horizon. 

4 Static and rotating equilibrium plasma configurations 

In this section we specialize to the study of stationary axisymmetric rigidly rotating plasma conflg- 
urations in a flat spacetime and discuss the corresponding phase diagrams of solutions. These also 
represent the phase diagrams for the dual black objects. We start with a general description of the 
system we want to study, and then we obtain the differential equation for the shape of equilibrium 
plasmas. We find the static plasma configurations in any dimension n, and the rotating plasma 
configurations for n = 1 and discuss the resulting phases. 

4.1 General formulation of the stationary and axisymmetric case 

The shape of stationary and axisymmetric plasma configurations is determined by the Young- 
Laplace equation. The background geometry is d = n-|-3 dimensional fiat spacetime Mt x x 

*Namely, p+P ~ p in the non-relativistic limit and the extrinsic curvature has time derivative terms that disappear 
in the non-relativistic limit. These time derivative terms are however vanishing when dealing with stationary solutions. 
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(with n > 1), and we choose to work in coordinates such that the metric reads 

ds"^ = -dt^ + dr^ + r^{d9^ + sin^ 6 dcp"^ + cos^ 6 dnl_2) + dz'^ , (4.1) 

where 6 G [0, vr/2], cj) £ [0, 2-k) and dQ'^_2 is the metric of the unit (n — 2)-sphere. For n = 1 ah the 
results can be reproduced by setting 6 = it/2. 

We consider rigidly rotating configurations with fluid velocity 

u^ = l {6^ + o;^<5%) , j={l-g^^o,l) , (4.2) 

with constant angular velocity a;^. Note that a boost uJz as well as rotation in other planes (if 
allowed by the dimension of the spacetime) could be easily added to our discussion. 

These plasma configurations satisfy the equilibrium equation of state (j2.20p . Their entropy S 
and conserved charges E, J are given, respectively, by (|3.3p and (j3.4p with the relevant Killing 
vectors associated to (14. ip being S, = dt and X(t> = d^. We will look for surfaces of revolution 
invariant under the action of these vectors. 



4.2 Static plasma lumps 

Our discussion will be very succint since these configurations were studied in [20] and we have 
rederived their main results. However, since we have an equation of state for the fluid we are able 
to discuss the properties of the solutions in a different manner: we display the phases in an entropy 
vs energy diagram. These quantities map directly into the entropy and energy of the dual black 
holes and black strings in SS-AdS. Therefore a direct comparison is possible to the diagrams that 
have been obtained for the black hole and black string phases in vacuum gravity. 

For a profile of the form r = R{z) the equations follow easily from the Young-Laplace equation 
()2.6p which, as we have discussed extensively, requires the mean curvature K of the surface of a 
static fluid to be constant. The equilibrium equation of state ()2.20p fixes this constant to 

A- (4.3) 

There are three families of solutions to this equation: 

• Uniform plasma tubes (UT): These solutions have constant radius r = Rq along the circle, so 
K = n/Ro. These are essentially configurations of a plasma ball ind — l = n-|-2 dimensions, 
times a straight line z, and have pressure and density 

P = n^, p = n(n + 3)^ + (n + 4)/)o. (4.4) 



Plasma balls (B): These solutions are characterized by R{z) = y^i?^ — , i.e., they describe 
a (n + l)-dimensional sphere of radius Rq, and extrinsic curvature K = {n + l)/Ro, so 

P={n+1)^, p=(n + l)(n + 3)-^ + (n + 4)po. (4.5) 
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Figure 1: Entropy versus energy phase diagrams. Tlie entropy is normalized to be one for uniform plasma 
tubes. Uniform plasma tubes, plasma balls and non-uniform plasma tubes correspond to squares, circles and 
triangles respectively. The left diagram, plotted for n ~ 1, is representative for n < 7, and shows that the 
preferred configuration is the plasma ball for small energy, and the uniform plasma tube for higher energy. 
The unstable plasma ball branch merges with the non-uniform plasma tube branch in the merger point. 
Above the critical dimension, for n > 8 (the right diagram corresponds to n = 12), the non-uniform plasma 
tubes become stable and become the preferred configuration for an intermediate range of energies. The left 
figure should be compared with Fig. 3 of [57] for the phases of non-uniform black strings. To plot these 
diagrams we used the values a — 1, po ^ and cr = 1. 

• Non-uniform plasma tubes (NUT): For n ^ 0, there is a third family of static equilibrium 
solutions that describe non- uniform plasma tubes. These configurations are found by solving 
numerically a first integral of the equation for the mean curvature. 

We analyzed these solutions numerically, and constructed the corresponding phase diagrams in 
the E — S plane. These are shown for n = 1 (left) and n = 12 (right) in Fig. [TJ The preferred 
configuration, at a given energy, is given by the one with highest entropy. There is a critical 
dimension = 8 for which the qualitative behavior changes. For n < 7, the entropically favored 
configuration is the plasma ball for small energy, which turns into a uniform plasma tube through 
a first order phase transition as the energy is increased. The non-uniform plasma tube branch 
has always a lower entropy (Fig. [H left). However, for n > 8, the uniform plasma tubes benefit 
entropically from generating non- uniformities in an intermediate range of energies (Fig. [H right). 

The case n = is special, since the fiuid lives in a two-dimensional space and the uniform plasma 
strips and plasma disks are the only solutions with constant mean curvature. Static plasma strip 
solutions have vanishing fiuid pressure: since the section of the boundary at constant z consists of 
just two points, there is no extrinsic curvature of the surface that must be balanced by pressure. 
This means that capillarity does not exert any force in the direction orthogonal to the plasma strip. 
This is consistent with the fact that the surface area does not change when the tube radius changes. 
In this case, the configuration is in a stable equilibrium, and no non-uniform plasma strip solution 
exists. 

In the next section, we shall obtain a better understanding of these features in terms of the 
mechanical stability properties of the uniform plasma tubes. Indeed, the existence of the non- 
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uniform plasma tube branch can be inferred from the fact that the uniform plasma tube changes 
its stability properties for a critical energy, at which it develops a static unstable mode. This 
critical unstable mode implies a bifurcation point to a new branch of static solutions that describe 
the non-uniform plasma tubes, and is associated to the Rayleigh-Plateau instability. For n = 0, 
the corresponding threshold mode has infinite wavelength, so there cannot be any non-uniform 
branches connected to the uniform strip. 

Note also that one can have an arbitrary number of coexisting disconnected plasma lumps. 
In the plasma fluid approximation, which is valid as long as they are far apart compared to the 
mean free path of the fluid, they do not feel each otheo In particular, other phases in which 
the fluid is distributed in multiple plasma balls can be present, but they have subleading entropy 
with respect to the single plasma ball in the phase diagram. Also, one can construct new branches 
of non-uniform plasma tubes, known as copy solutions, by unwrapping the original non-uniform 
plasma tube along the compact z direction, and changing its periodicity. These new branches can 
be easily obtained through scaling arguments [28]. Since they are subleading in entropy we will not 
consider them. 



4.3 Rotating plasma lumps 

The problem of finding general equilibrium solutions in dimension n is considerably complicated 
by the presence of rotation, which breaks the SO{n + 1) symmetry of the sections at constant z 
to SO{2) X SO{n — 1) and introduces a dependence on the angle in the profile of the lumps. 
Then one is forced to solve a partial differential equation. However, this problem is not present 
when n = 1, where the transverse spheres are actually circles so there is the same symmetry as in 
the static case and the equations reduce to ODEs. Therefore we restrict our analysis of spinning 
lumps to n = 1. Even if there may be phenomena of higher-dimensional tubes that we miss by 
this restriction, the case of n = 1 still exhibits qualitatively new dynamics relative to the static 
situation. 

We represent the boundary of the spinning lumps in terms of a height function /i(r)j£| 

f{r,z)=z-h{r) = (4.6) 

Then, use of the equilibrium equation of state (I2.20p for the pressure and of the Young-Laplace 
equation ()2.6p yields the ODE for a stationary axisymmetric plasma configuration, 

d f \ , P* 2 2N-5/2 Po „ 

— —== + —r [I — r LOj, r = 0. (4.7) 

Some solutions to these equations were already studied in detail in [T31 E], namely, balls, 
pinched balls and rings. Now they are constrained to fit along the z direction, which we are taking 
to be compact. Two other solutions, the rotating uniform tube and the uniform 'hollow tube' can be 

®For the dual black holes, this means separations larger than the cosmological length scale l/V— A, and their 
mutual gravitational interaction would involve massive Kaluza-Klein modes in the vacuum phase, which is subleading. 

^For the uniform tubes, the simplest boundary parametrization is instead /(r) — r — Ro = and the Young-Laplace 
equation gives pc as a function of po, Ro-, oj^. 
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readily constructed out of them by simply takina; the rotating plasma disk and the plasma annulus 
that appear in one less dimension and trLslating them uniformly along J To make contact with 
this previous work we follow partially their notation and define the dimensionless variables, 

uj^ = ^^, v = uj^r, H{v)=ujMr), k = \—, L = —L , (4.8) 

where L will be used to represent the length of the uniform tube with radius Ro and L the corre- 
sponding dimensionless quantity. Note also that reality of 7 = (1 — requires < v < 1. 
Integrating twice (14. 7p we get the function H{v) that describes the profile of the spinning lumps, 

H{v)= I '^^^^ jdx; with f(v) = 2k-?>{v'^ + 2c)-i-^, g{v) = Qui^v-f'^ . 

Jvo {g{x) — f{x) ) 2 

where c is an integration constant. Plasma rings have also a inner surface: its profile is described by 
a similar expression obtained by multiplying the first relation in (j4.9p by —1 and replacing Vq — > v\ 
(with v\ being the inner velocity; for details see [H]). 

To make contact again with [T3] we define the dimensionless energy, angular momentum and 
entropy as 

^ = |i4, 7=^, S = ^%. (4.10) 

47ro-^ 207rcr4 207ra 5 cr^ 



where E, J, S follow from (j3.4p and (j3.3|) . It is also useful to introduce the following functions: 

Etivo) = -^{2k-f',{vl + 2-2%'')+3vo{vo + ^^)) 

Mvo) = ^{2^-'-2 + 5vl) , 5,(„,) = _^(^3_i)_ (411^ 



and 



Ebia, b) = ^ [ dv (vH{v) {k{A + v'^W + l) + ^4,v^/lVW(vf 
~ k ~ k^l^ 

Jb{a,h) = ^ v^j'^H{v)dv, Sb{a,b) = v-f^H{v)dv. (4.12) 

J a J a 

We now discuss the seven families of axisymmetric spinning plasma lumps that we found to be 
solutions of ()4.7p (we present the ones with non-trivial profiles in Figs. [21 [3l and [H The figures 
should be rotated along the vertical axis v = 0): 

• Uniform plasma tubes (UT): These have constant radius Ro along their length L. They 
have u)A, = Vo{kj^ — 1), where 7^ = 7! . Their energy, angular momentum and entropy are 

E = Et{vo), J = Jtiyo) and S = St{vo) as defined in ()4.1ip . 



''Observe that this is not possible for pinched balls since they do not exist for n = 0. 
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Figure 2: Figures of equilibrium of a rotating plasma tube. From left to right, an ordinary non-uniform 
rotating plasma tube {vm = -5, Wo — -8, k — .628024), a pinched plasma tube = .201016, Vo = .918, 
k — .15) and a plasma tube on the verge of splitting off a plasma ring {vm — .201016, Vo — .9133, k — .15). 
The parameters have been chosen such that the three tubes share the same periodicity. 



Uniform hollow tubes (UHT): These have two constant radii, outer Rq and inner R\, along 
their length L delimiting the tubular region inside which the fluid lies. The Young-Laplace 
equation on the inner and outer boundaries fixes a constraint between these two radii, which 
can be written as oja, = Vo(k^o — 1) = fi(l — k^f), where 7oi = 7! • The energy, angular 

momentum and entropy are E = Et{vo) — Et{vi), J = Jt{vo) — Jt{y\) and S = Stivo) — St{vi), 
with the definitions in (j4.1ip . 

Non-uniform plasma tubes (NUT): These tubes have a radius that varies from a non-vanishing 
minimum, Vm, tip to a maximum, Vo, as we go along their length. At both these extreme 
points we have H'{v) = —00, and in between them H'{v) is always negative (so the point in 
between the two extreme points where H"{v) changes sign has H'{v) < 0); see Fig. [2la. In 
the critical NUT configuration, H'{v) = at the infiection point and the NUT and pinched 
NUT (discussed next) meet in what we call the critical NUT. These boundary conditions 
imply the following relations between the parameters, 

^ /1 2 ^-3/2 1 2 



[^k[{l-vly"'-{l-viy"')-^vl-vl)\. (4.13) 



1 v„ 
6wo 



The energy, angular momentum and entropy of the NUTs are given, respectively, by -E = 
Et{vm) +Eb{vm,Vo), J = Jt{vm) + Jb{vm,Vo) and S = St{vm) + Eb{vm,Vo) where we use the 
functions defined in ()4.1ip and ()4.12p . 
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Figure 3: Figures of equilibrium of a rotating plasma ball. From left to right, an ordinary rotating drop 
{vo — .7, k = 2), a marginal rotating drop (vq = .1, k = 21), a pinched plasma ball {vo = .5, fc = .83) and a 
pinched ball on the verge of splitting of a plasma ring [vq = .7965, k = .5). 



0.2- 




Figure 4: Figures of equilibrium of a rotating plasma ring. On the left, a thin plasma ring (wj = .2, Vo = .6, 
fc = 9) and on the right a fat plasma ring (wj = .1, = .4, fc = 2). 

• Pinched non-uniform plasma tubes (pNUT): These tubes also have a radius that varies from 
a non- vanishing minimum, Vm, up to a maximum, Vq- At both these extreme points we have 
H'[v) = —oo. But, contrary to the NUT solutions, H'{v) now changes sign twice as v runs 
from Vm to Vo (so the point in between the two extreme points where H"[v) changes sign has 
now H'{v) > 0), see Fig. Ob. These boundary conditions still imply the relations (14.131) and 
the conserved charges of the pinched NUTs are given by the same functional defined for the 
NUT in the previous item. As the rotation increases we reach a configuration where we have 
a pNUT and a plasma tube on the verge of splitting off a plasma ring; see Fig. [2lc. 

• Plasma balls (B): These plasma configurations are oblate spheres. At the extreme points 
they satisfy H'{v = 0) = and H'{v = Vq) = — oo. Moreover, H'(v) is always negative at 
the boundary as we go from f = up to u = Uq; see Fig. Oa. The condition H'{0) = fixes 
c = k/3 and H'{vo) = — oo requires f{vo) = g{vo) which fixes tJ^ = (2fc7o — 2k — Sv^)- 
Finally, the condition H'(v) < in the interval (0, fo) implies k > 1. For k = 1 the critical 
ball configuration is reached (see Fig. Ob). Here, the plasma ball joins the pinched ball branch 
discussed next. The conserved charges are given, respectively, hy E = Eh{0, Vq), J = Jb{0, Vq) 
and S = Eb{0,Vo) where we use the functions defined in (j4.12p . The balls must have a 
distance between their poles smaller than the compact dimension length. The same applies 
to the solutions we discuss next. 
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• Pinched plasma balls (pB): These balls also have H'{0) = and H'{vo) = — oo which again 
requires c = k/3 and tD^ = g^- (2/07^ — 2k — 3Wq). But contrary to the plasma balls, H'{v) 
changes from positive to negative sign in the interval (0, fo); see Fig. [3jc. This requires 
< A; < 1. However not all k in this interval give a physical solution. They must further 
satisfy H[v = 0) > 0. At the critical situation where H{{)) = the pinched plasma ball 
branch joins the fat plasma ring phase discussed next (see Fig. Od). The conserved charges 
of the pinched balls are given by the same expressions as for the balls. 

• Plasma rings (R): These annulus configurations are characterized by H'[v^) = +00 and 

o) — —00 at the inner and outer boundary, respectively, and for having H[v.J — H{vq) — 
0. This requires that /(vj = —g{v-) and f{vo) = g{vo)- For fixed outer boundary velocity, 
Vo = uj(j,Ro, there are two possible inner velocities satisfying these conditions. The smaller 
one describes the fat plasma ring (that joins the pinched ball branch at the critical configura- 
tion where v. 0) and the larger v. describes the thin ring (see Fig. The two rings meet 
at a regular critical ring when the two aforementioned v.^s are equal. The energy, angular 
momentum and entropy of the rings are given, respectively, by -E = Eb{v^,Vo), J = Jb{v^,Vo) 
and S = Eb{v.,Vo) with functions defined in (j4.12p . 

There also appears the possibility of non-uniform hollow tubes (NUHT), i.e., hollow tubes with 
an inhomogeneous profile along z, but we have not obtained direct evidence for them. If they 
exist, it would be interesting to increase their rotation and investigate their possible pinches and 
connections to other phases. 

Having identified these families of plasma solutions, it should be interesting to represent their 
properties in a phase diagram where we fix the periodic direction length and the energy of the 
solution, and let the entropy vary with the angular momentum (alternatively, we could also fix L 
and J and represent the solutions in the S{E) diagram). This requires substantial and delicate 
numerical work since this is a shooting problem on three parameters. The curves for uniform tubes 
and balls with a fixed conserved charge can be obtained relatively easily but the numerical search 
of the non-uniform tubes is considerably more difficult. We leave this for future work. 

5 Rayleigh-Plateau instability of plasma tubes 

Non-relativistic fluid jets and liquid bridges in between two disk plates are unstable against the 
so-called Rayleigh-Plateau instability if their length is larger than their transverse perimeter. This 
instability is active even in the absence of gravity and surface tension is the crucial mechanism 
responsible for it: in the simplest case, capillarity forces conspire to reduce the surface area of the 
fluid. This is a long wavelength instability that makes the fluid cylinder pinch-off. 

In this section we show that uniform plasma tubes are naturally afflicted by the relativistic 
Rayleigh-Plateau instability and we will study in detail its properties. We start by getting a good 
understanding of the mechanism responsible for the instability noting that there are perturbation 
modes that decrease the plasma tube potential energy while keeping its volume fixed (so this 
statement is independent of whether or not the fluid is compressible, i.e., of the equation of state). 
Equivalently, in view of the results of section [Sj these modes are unstable because they increase the 
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entropy of the system while keeping its conserved charges fixed. In subsection [STTl this approach will 
allow to find the minimum length below which the uniform tube can be stable. Then in subsection 
15.21 we will find the dispersion relation for this instability. 

In the dual gravitational system, uniform black strings are expected to suffer from the grav- 
itational Gregory-Laflamme instability. It is also a long wavelength instability that makes black 
strings longer than the horizon radius pinch. We will see that the plasma surface tension and 
gravitational instabilities have similar properties. 

5.1 Unstable threshold mode from entropy or potential energy variation 
5.1.1 Threshold mode for static and boosted plasmas 

Consider an axisymmetric perturbation on a uniform plasma tube, 

r{z) = Ro + eRi cos{kz) + e^i?2 , e < 1 . (5.1) 

We take the tube length to be given by the Rayleigh-Plateau wavelength because that is where it 
might break, L = X = 27r//c. The volume of the perturbed tube, up to a second order expansion, 
is then 

_(n+l)/2 / „ _|_ 1 /„ X \ 



We require that the volume of the tube is kept fixed up to second order in the perturbation 
(equivalently, its conserved charges must be kept fixed) so this constrains R2 as a function of the 
lower order radius, 

= (5.3) 

The difference between the perturbed and unperturbed surface tension potential energy is 

2„("+i)/2 

where we used ()5.3p . Unstable modes are those that decrease the potential energy, i.e., surface 
area, of the fluid tube. So the condition Af/o- ^ for flxed volume requires 

kRoi^^/^, (5.5) 

For static plasmas the threshold unstable mode is then k = y/n/Ro and it increases with the spatial 
dimension. Finally, a boost along the z-direction has only a kinematical effect on this instability, 
Lorentz-contracting the threshold wavelength [29]. 



5.1.2 Threshold mode for rotating plasma and non-axisymmetric modes 

We now generalize the previous study to rotating and boosted uniform plasma tubes, as well as 
non-axisymmetric perturbations on it. Again, we will restrict to the n = 1 dimensional case, since 
for higher n, the profile of the drop has a polar angle dependence that prevents an analytical study. 
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The general perturbation on the uniform plasma tube reads, 

r{z, <p) = Ro + €Ricos{kz) cos{m(l)) + e'^R2 , e <C 1 . (5-6) 

This perturbation includes axisymmetric modes (m = 0) as well as non-axisymmetric ones (integer 
m 7^ 0). Again, we take the tube length to be equal to the instability wavelength, L = X = 27r/k. 
The volume of the perturbed tube, up to a second order expansion, is 

V = XttRI + ^ (vr {Rj + 8R0R2) + Rl ^^^^1^) + 0{e') . (5.7) 

We require that the volume of the tube is kept fixed up to second order in the perturbation 
(equivalently, its conserved charges must be kept fixed) so this constrains R2 as a function of the 
lower order radius, 

i?2 = -^, form = 0; R^ = -^, for m ^ . (5.8) 

The potential energy of the tube is the sum of the surface tension and centrifugal contributions 
as defined in section [3l Using the constraint (|5.8p , the total potential energy difference between the 
perturbed and unperturbed configurations At/ = AC/^ + AC/cf is then given by {5m = 1 if m = 
and vanishes otherwise) 

AC/ = ^ (1 + 5m) (y{k'Rl - 1 + m^) - (1 - Rlu:l - ul)-'''"^ R\^ + 0{e') . (5.9) 

The Rayleigh-Plateau instability is active for wavenumbers that decrease the potential energy of 
the fluid tube. So the condition AU ^ for fixed volume requires 

k'Rl + m'^l + (1 - Rl4 - u^l)-"' , (5.10) 

The inverse of the ratio p^uj'^R'^/a, that measures the competition between the surface tension 
and centrifugal effects, is often called the rotating Bond number or Hocking parameter. Rotation 
increases the critical wavenumber making the plasma tube unstable for a wider range of wave- 
lengths. Only axisymmetric modes are unstable in the static case but when rotation is added, 
non-axisymmetric modes can become also unstable for sufficiently high velocity. We will discuss 
the dual gravitational interpretation of these unstable modes in section [7l 

5.2 Rayleigh-Plateau dispersion relation 

In this subsection we want to address the stability of a uniform plasma tube when we perturb 
it. The dynamics of the perturbations is dictated by the hydrodynamic equations, subject to 
appropriate boundary conditions. 

Perturbations take the plasma away from thermal equilibrium and therefore viscosity and diffu- 
sion effects start to contribute. The energy-momentum tensor of the fluid includes now not only the 
perfect fluid and the boundary surface tension terms (j2.ip . but also a dissipative contribution. As 
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we will show, the uniform plasma tube is afflicted by the Rayleigh-Plateau instability and surface 
tension is the mechanism responsible for it. Viscosity and diffusion play no role on the activation 
of the instability and have a subleading effect on the dispersion relation, simply correcting the 
threshold wavelength and the time-scale of the instability. Hence, in our analysis we shall neglect 
the dissipation contribution to the fluid energy-momentum tensor, and simply comment at the end 
the subleading effects it introduces. 

To study the stability of rigidly rotating uniform plasma tubes, we consider a generic un- 
perturbed uniform tube with velocity (j4.2p . The precise expression for the unperturbed plasma 
pressure and density is fixed once we specify the equation of state. In particular, for the plasma we 
are interested in, these are specified by equations (12.160 and ()2.20p . The particular choice of the 
equation of state will however not be fundamental in our analysis. 

A generic perturbation on the uniform tube is described as 

= nfo) + (5n^ , P = P(o) + <5P , p = p^^-^+5p, (5.11) 

where we denote an unperturbed quantity by (5(o) , and the perturbation as 5Q. The perturbed state 
(15. lip must satisfy the relativistic continuity and Navier-Stokes equations, 1)2.31) and (12. 4p . The 
Young-Laplace equation (j2.6p provides then a boundary condition for the perturbed problem. After 
eliminating the 0**^ order terms using the unperturbed hydrodynamic equations, the continuity and 
the Navier-Stokes equations yield, up to first order in the perturbation, 

u^^^^d^dp + 5n'^a^P(o) + (/0(o) + ^'{0))V^<5n'^ + {6p + 5P)V ^u^^^^ = , (5.12) 

(P(0) + ^(0)) (-^n^V^n^o) + nj'o^V^^n'^) + {5p + 5P^^^V ^u\^-^ 

+ (g^^ + nfo^tx^o)) V^5P + (^n'^nj,) + u^^^Su^) V^P(o) = . (5.13) 

The density and pressure perturbations are not independent. They are related by the equation of 
state (|2.16p . valid also out of equilibrium. In particular, perturbation of the first relation in (j2.16p 
yields 

5p={n + 3)6P . (5.14) 

Since we can expand any perturbation in a Fourier series, we restrict the analysis to a generic 
mode and consider a perturbation that disturbs the boundary surface according to 

r = R{t,z,^), R{t,z,^) = Ro + ee'^^e'^^^'"'^ , e ^ Ro , (5.15) 

where Rq is the unperturbed radius of the uniform tube (note that for the rotating n > 1 case 
it is ^-dependent). Positive uj describes an instability with wavenumber k. These modes break 
axisymmetry if m 7^ 0. The unit normal of ()5.15p is 

n, = \n\-'{-R[5;, + 5;;-R'^5^-R'J^) , \n\ = h - R[^ + ^ + R'M . (5.16) 

Naturally, we look for perturbations of the fluid quantities that have the same form as the boundary 
disturbance, 

6Q{t, r, z, 0) = 5Q(r)e'^*e*^^+^'"'^ , 6Q = {6u'', 6P, 6p} . (5.17) 
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The perturbed hydro dynamic equations ()5.12p and ()5.13p must be supplemented by appropriate 
boundary conditions. The first one demands normal stress balance on the boundary. This means 
that the pressure perturbation that solves (|5.12p and (j5.13p must also satisfy the third perturbed 
hydrodynamic equation, namely the equation that follows from perturbing the Young-Laplace 
equation (j2.6p . 



BC I : 6Pl, ~ a 

I bdry 



K\ — K\ 



< > J \R{t,z,4>) \ < > J \Ro 



, (5.18) 



where on the rhs we evaluate the expression at the perturbed boundary r = R{t, z, (jj) defined in 
(|5.15p and subtract the unperturbed contribution evaluated at r = Ro- For our plasma, — P^^ 
is obtained from (j2.20p . and the extrinsic curvature is given from its definition (|2.6p using the unit 
normal (j5.16p . 

The second boundary condition is a kinematic condition requiring that the normal component 
of the fluid velocity on the boundary satisfies the perturbed version of (12. 2p . u'^^^^ 6n^ + 6u^n^/^^ = 0, 

where 6n^ = n^\j^^^ z <p)~ ^^^^ ^'^'^ unperturbed normal is n}?^ = n^|^ = This ensures that 
the velocity perturbation leaves the fluid confined inside the boundary. This boundary condition 
then reads 



BC II : ^^^Ibdry - (1 - ^I'^l - ^ ('^ + ^"^'^'^ + ^^^-) ee'"*e^^"+^'"<^ . (5.19) 



We have now all the ingredients needed to find perturbations of the form (j5.15p and (|5.17p 
that might develop an instability on the tube. We will find that there is indeed a long wavelength 
instability known in non-relativistic systems as the Rayleigh-Plateau instability (see for example 
[30]). It afflicts our plasma tube when, roughly, its length is larger than its transverse radius. 
We will first analyze this instability and its dispersion relation uj{k) for a static uniform tube in 
any dimension n. This is the simplest case where the crucial ingredients necessary to activate the 
instability are present. Then we will study the changes introduced in the instability properties by 
rotation and boost, as well as the effects of non-axisymmetric modes. This will be done in the n = 1 
case, where the analysis can be done analytically, but we expect the results to hold qualitatively 
also for n> \. Finally, we will discuss briefly the subleading effects that viscosity would introduce 
in the dispersion relation. 



5.2.1 Perturbations of static tubes in any dimension 

Consider a static uniform plasma tube with unperturbed radius Ro and u^'q^ = (5^^, in [n + 3) 
flat spacetime dimensions, and deform its boundary by axisymmetric perturbations {m = 0) of the 
form (|5.15p . The perturbed continuity and Navier-Stokes equations, (|5.12p and (|5.13p have unstable 
modes that correspond to the Rayleigh-Plateau instability, as we show in the sequel. Solving (j5.13p 
with perturbations (I5.17p . we find that the non- vanishing 5u^ are (at higher order other components 
will arise) 

5u^{r) = -u;-i (p(o) + P(o))"' , 5u^{r) = - (p(o) + ^'(o))"' 5P{r) , (5.20) 
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which, replaced in the continuity equation ()5.12p . give 



1 



-p^dP{r) =0, p = kil + {n + 3)-^] . (5.21) 



dr"^ r dr ' \ k"^ ^ 

The uj term in the definition of p arises from the first term in (|5.12p proportional to 5p and use 
of the equation of state ()5.14p . Equation (|5.2ip is a modified Bessel equation, whose solutions are 
the modified Bessel functions of the first kind I+n—l{pf) a-nd second kind Kn^ipr). Of these, 

2 2 

Tl — 1 

I_n^(pr) and K n-i (pr) diverge as r ~ as r — > and we discard them. Therefore, the regular 

2 2 

solution of (j5.2ip at the origin is 



n-l 



5P{r) = Ar 2 /^-i (yr) , 



2 



5vl'{r) = „ X r "2' fp !'„_-, (pr) - ^^^^^ In^{pr) 

w(P(0) +-P(0)) V ~ 2r 2 

6u'{r) = rr-V/^(pr), (5.22) 

^ (P(o) + ^(0) ) 2 

where we defined Il{x) = dxlu{x), and where A is a constant that is fixed by the boundary condition 
(|5.18p to the leading-order value 



n — 1 2 

^ - e^T^^V^ (^^^o + ^^Rl - n) ■ (5.23) 
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With this knowledge, and using the relations /^(x) = Iu+i{x) + ^I^{x) and (j2.20p . the boundary 
condition (j5.19p yields the desired dispersion relation oj = uj{k), which reads 



2 

A plot of this dispersion relation a; (A;) for several values of the dimension n is shown in Fig. [5l To 
obtain these plots we have needed to specify the value of ^ . This amounts to considering plasma 
tubes of a specific width Rq relative to the length scale cr/p*, which is essentially the mean free 
path of the plasma. In principle we should vary this width and study the dispersion relation as a 
function of it. The shape of the dispersion relation is indeed modified as ^^"^ grows to values of 
order one. However, the validity of the fiuid description requires that we consider small ^ (see 
section [6|). It turns out that in this case the qualitative aspects of the dispersion relation hardly 
change for different values of ^ and the shapes shown in Fig. [5] are generic. For definiteness, we 
have taken 

^^10-^ (5.25) 



/9*i?o n + 6 

so as to produce the plot. But we insist that the dimension-dependence in this choice, and hence 
in the curves in the plot, is completely arbitrary: it only amounts to choosing tubes of a particular 
width in each dimension. 
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Figure 5: Plot of the dimensionless dispersion relation ui{k) for the Rayleigh-Plateau instability in a static 
uniform tube for several spacetime dimensions d = n + 3. The instability strength and threshold wavenumber 
increase as the spacetime dimension grows. Intuitively, this is because the "same" perturbation decreases 
more the surface area of the plasma boundary (increases more the plasma entropy) as the dimension increases. 



Note also that the most unstable mode satisfies loRq 



max Pq Ro 



and is thus within the regime of validity 



of the hydrodynamic description (see section [Q. This plot is to be qualitatively compared with Fig. 3 of [S] 
for the GL instability. 



The Rayleigh-Plateau instability is active when w > 0. The instability gets stronger as the 
spacetime dimension increases since the threshold wavenumber, as well as the wavenumber and 
frequency where the maximum instability occurs, increase when n grows. 

There are two special points that we can read easily from (|5.24p . One is the point where 
the curves cross w = 0, which gives the threshold wavenumber below which the Rayleigh-Plateau 
instability is active: kRo = \/n. So unstable modes satisfy condition (j5.5p . confirming the threshold 
mode that we obtained previously from entropical or energetical arguments in subsection 15.11 The 
second point is at A; = 0, and one can expand the Bessel functions in (|5.24|) to check that uj at 
A: = always vanishes, independently of -jp^- 

5.2.2 Rotating and boosted plasma lumps 

We want to consider now how rotation and boost of the fluid modify the dispersion relation, and 
study the effect of non-axisymmetric modes (m ^ 0). We shall restrict to the n = 1 case where an 
analytical approach is viable. Take a uniform plasma tube rigidly rotating with velocity (j4.2p and 
with unperturbed constant radius Rq- The unperturbed state satisfies, once the equation of state 
for the fluid is introduced, equations (|2.16p and (|2.20p with n = 1. 

Disturb now the plasma boundary with perturbations of the form (jS.lSp . To get the velocity 
perturbations we have to solve (j5.13p assuming perturbations of the form (j5.17p . One finds, after 
use of (|5.14p . (in (|5.26p - (|5.28p we restore the velocity of light factor c to make the perturbative 
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expansion that we do next more clear) 

„2 



5u\r) 

Su'l'ir) = ^ 




37 r ■ 



5P{r) 



5p* ^2 _^ 4^2^ 
4 k 

6P{r) 



6P{r) 



47 

5/3* 



ZI7 



where we defined 



w = u + imtOfp + ikuJz 



X = 1 + 



3a;^72r 



6P{r) 



(5.26) 



(5.27) 



Note that (j5.26p satisfies the perturbed version of u^u^ 
(|5.26p can now be replaced in the continuity equation, 

d6u 



c^, namely u^^^6u^ = 0. Equations 



, +-5u'' + im 5u'^ + ik 5u'' + w (5n* + 47^r ^5u'' + ^^^5P = , 
dr r 5p*c^ 



(5.28) 



r 

As will become clear, it is not possible to get an analytical expression for the dispersion relation 
unless we work in the small velocity regime, uj^pRo, uJz ^ c. Since the Rayleigh-Plateau instability is 
already present in the static case and we just want to find if adding velocity to the solution increases 
or decreases the instability strength and threshold unstable mode, we will solve the hydrodynamic 
equations in this regime (we will further represent higher order terms of the velocity generically by 
0{uJi)). This restriction is also justified by the fact that the regime of validity of the hydrodynamic 
description and thus of our results is restricted to small rotation rates, as we will argue in section 
[H So, to leading order (j5.26p reduces to 



6u\r) 
<5u"(r) 



0(lo 



5/j* vu^ 



4a; 



2 I Or oP{r) + ——dP{r^ 



w 



1 



'^'^''dr 5P{r) - '-^5P{r) ] + Oii^'f) , 



5p* tu^ + 4a;2 r \ w 

5u''{r) = 

The perturbed continuity equation yields then to leading order in the velocitylf] 

4a;2\i 



d'^5P{r) ^ 1 d5P{r) 



dr^ 



dr 



m 



<5P(r) + 0(w/) = 0, r] = k[l + 



1 + 



(5.29) 



1/2 



(5.30) 



^Notice that the last term in (|5.28p has a contribution ijjjf? that we keep because our expansion is only on the 
velocities but not on the frequency. This contribution is responsible for the u)^ /k^ term in -q. 
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This is again a modified Bessel equation with solutions Im{w) ^^^d Kmirir) (note that if we had 
kept higher order terms we would not be able to solve analytically the differential equation). We 
discard the Kmirjr) solution that diverges at the origin and thus obtain the regular solution 



6P{r) ~ Al^niw) , 



+ 4a;j 



(5.31) 



The final expression for du'l', (that we do not need) can then be obtained from (|5.29p . 
The boundary condition (|5.18|) fixes the ratio A/e to leading order as 



A 



ctR- 



IfnivRo) 



k^Rl + oj'^Rl - 1 + 



5 P*^lRl 



1 — uAR^, 



UJ. 



a 



-7/2 



(5.32) 



Use of this ratio on the boundary condition (I5.19P yields finally the leading dispersion relation uj{k) 




ImillRo) 



+ m [ 1 



2o3 



W J 

,2\-7/2 



2d2 



(5.33) 



A plot of this dispersion relation uj{k) for different values of w,^ (and = m = 0) is shown in 
Fig. [6l The dispersion relation (j5.33p . valid in the small velocity regime, is expect to be a good 
approximation up to velocities of the order Vq ~ 0.1 — 0.5. The Rayleigh-Plateau instability is 
active when uj > 0. The marginal unstable mode with uj = satisfies condition (|5.10p . confirming 
the threshold mode obtained in subsection 15.1.21 As the rotation is increased, the threshold 
wavenumber, as well as the wavenumber and frequency where the maximum instability occurs, and 
the instability grows stronger. Fixing a;^ ^ 0, increasing ujz also makes the instability stronger. 
Finally note that contrary to the static case where only axisymmetric modes were unstable, now 
non-axisymmetric modes can also become unstable for sufficiently high rotation. 

Note that we decided to keep the higher order terms in velocities in the boundary condition 
(|5.32p and in (j5.33p instead of doing a Taylor expansion. The reason being that from the en- 
tropic/energetic computation we get the threshold condition ()5.10p in the marginal unstable case. 
The computation leading to this condition is exact to all orders in the velocity, so we do not expand 
()5.32p that is responsible for the unstable wavenumber cut-off in (j5.33p . 

In section [6] we will discuss the regime of validity of the hydrodynamic description. Our results 
are valid for {ujRo,kRo} ^ J^T' There is a wide range of parameters for which this condition is 
satisfied. In particular, the physically relevant most unstable mode of Figs [5] and [6] fits very well 
in this regime. Near the extreme points of the dispersion relation where a; = the results must be 
read with caution. 



5.2.3 Viscosity 

We have neglected in our analysis the effects of viscosity, which is absent from stationary solutions 
but should play a role in time-dependent processes such as the Rayleigh-Plateau instability. 
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Figure 6: Plot of the dimensionless dispersion relation ui{k) for the Rayleigh-Plateau instability in a rotating 
uniform tube for dimensionless angular velocity = 0; 5 x 10"''; 10"'^, where Vo = oJ^Ro- The instability 
strength and threshold wavenumber increase as the rotation grows. Intuitively, this is because the centrifugal 
force scales with the radial distance and thus it is bigger in a crest than in a trsough of the instability. The 
numerical data corresponds to take 5-^^ — 10^^, = and m = 0. The dispersion relation (|5.33p . valid 
in the small velocity regime, is expect to be a good approximation up to velocities of the order Wo ~ 0.1 — 0.5. 
Note also that the most unstable mode satisfies ujRo I ^ — %- , and is thus within the regime of validity 

" I max po -Ro ' " 

of the hydrodynamic description (see section [6]). This plot should be compared with Fig. 1 of 40J for the 
GL instability of rotating black strings. 

These effects are well understood in the the non-relativistic Rayleigh-Plateau instability, which 
has been the object of exhaustive theoretical, numerical and experimental study. It does not seem 
unreasonable to assume that the effect of viscosity in the relativistic and non-relativistic instabilities 
is qualitatively similar. 

Viscosity plays a subleading role in the activation of the RP instability, although the instability 
can get considerably weaker if the viscosity is very high. Generically, it increases the wavelength 
of the most unstable mode, and weakens the strength of the instability [30]. Viscosity has been 
experimentally seen to play an important role at later stages in the time evolution of the instability, 
when the fluid tube pinches off, as we discuss in sec. [7] 

Dissipation increases considerably the technical challenge of solving numerically the hydrody- 
namic perturbed equations. This is because without it, these are ODEs with boundary conditions, 
while when dissipation is present we must deal with a coupled system of second order differential 
equations for the pressure and velocity perturbations, which cannot be studied analytically and 
requires advanced numerical treatment. 

6 Regime of validity 

We now discuss the limits of validity of the different approximations we have made in our derivations 
in the previous sections. 

First of all, for the fluid description of the deconfined plasma to be accurate the thermodynamic 
quantities must vary slowly over the mean free path i^fp of the fluid, which is of the order of the 
mass gap of the theory, or equivalently of the order of the deconfinement temperature. That is, all 
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length scales A in the fluid must be 

A » T^^ ~ — . (6.34) 

For the same reason, thermal fluctuations must remain small. A good estimate for the valid regime 
is obtained when the maximum fractional rate of change of the fluid local temperature is small 

5T, 



To 



^rlHTL. <<1. (6.35) 



I max ' ' I max 



This occurs for 

- — ^ «; . (6.36 

For plasma balls and rings this condition is satisfied for a wide range of lumps as long as we are 
away from extremality where the temperature vanishes [13]. For plasma tubes the angular velocity 
iOcf) increases with Ro so the above condition is fulfilled as long as the tubes do not rotate too fast. 
This does not affect our analysis in section 15.2.21 since, in order to have analytic control on the 
equations, we only considered the regime of small rotation. This condition also imposes restrictions 
on the validity of the Rayleigh-Plateau instability analysis. To guarantee that the thermodynamic 
quantities of the fiuid vary on a length scale small when compared with the mean free path we 
must demand that the Rayleigh-Plateau unstable frequencies and wavenumbers satisfy 

{ujRo, kRo} » . (6.37) 

Since the most unstable mode dominates the instability we must guarantee that this condition is 
satisfied in the vicinity of the maximum of the dispersion relation. We find that for ^^'^ < 10^^ 
this condition is satisfied, and things get better as becomes smaller. Recall that for the plots 
of Figs. [5] and [6] we took ^^"^ ~ 10~^, well within the validity range. 

Second, the interface between the confined and deconfined phases, i.e., the fiuid surface, has a 
finite thickness of the order 1/Tc [E], and therefore the delta- like surface approximation we used is 
valid provided that the curvature of the surface is small with respect to l/T^. So, for plasma balls 
and rings the analysis is valid when the boundary radii i?o! Ri, and Rq — Ri are large compared to 
TjT^. This is the case if the plasma energy is large and away from the extremal configurations [13j . 
For plasma tubes the requirements are, analogously, Ro, Rm, and Ro — Rm much larger than T~^. 

Finally, we neglected the dependence of the surface tension on the temperature and other 
thermodynamic quantities of the fiuid. For consistency we must then demand that on the boundary 
between the confined and deconfined phases the temperature of the plasma must remain everywhere 
close to the critical temperature Tc. The values of T and Tc at the boundary surfaces remain close 
for a large range of energies and angular momenta (see Fig. 9 of [l3]), both for balls and rings as 
well as for plasma tubes, always away from the extremal limits. 

7 Rayleigh-Plateau and Gregory-Laflamme 

The dynamics of lumps of fiuid that we have analyzed should admit, via the reasoning in [7], 
a mapping to the dynamics of black holes near the infrared bottom in SS-AdS, in the limit in 
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which the black hole size is much larger than the thickness of the domain wall that connects it 
to the confining 'AdS soliton' vacuum. In this case, the correspondence between black holes and 
lumps of fluid can be established and the solutions of the fluid equations can be regarded as actual 
gravitational calculations in a controlled perturbative (in boundary derivatives) expansion to the 
Einstein equations in AdS, reformulated in an appropriate way. 

This implies that our results for the dynamics of plasma tubes in d dimensions can be regarded 
as describing gravitational physics of black strings in SS-AdSrf+2 spacetime. However, such black 
string solutions have not been constructed directly yet and therefore a comparison cannot be made 
between results obtained using dual-fluid techniques and more conventional gravitational analyses. 

Nevertheless, it is remarkable that the dynamics of fluid tubes from our analysis is strikingly 
similar to what is known for black strings in asymptotically flat space, for which a fairly large 
amount of results are available. It has been pointed out that all the different plasma ball, pinched 
ball and plasma ring phases found in [13^ [T^ in d dimensions have black hole, pinched black hole, 
and black ring analogues in asymptotically flat d + 2-dimensional gravity [311 |32l [33] • 

This parallelism extends to most aspects of the solutions we have discussed. Ref. [20] already 
observed that the phase diagrams for static plasma lumps (shown here in Fig. [T|) resemble the 
corresponding phase diagrams for black strings and black holes localized in a Kaluza-Klein circle 
[Ml dS]. Let us add that if we compare our Fig. [Ha with Fig. 3 of [2^, the two figures can be 
strikingly superposed. Note that in Fig. 3 of [27J the gravitational merger between the non-uniform 
black string and localized black hole branches is not shown because the available numerical code 
breaks in this region [3l]. It is however conjectured that the two branches do indeed converge into 
a topology-changing merger with a conical structure [35]) ^-i^d there is numerical evidence for this 
[36]. The approximations for the fluid description become invalid close to this transition point (as 
do, too, the gravitational ones), but certainly a topology-changing transition seems to take place 
when this is neglected and the calculations are pushed until pinch-off. The similarities between the 
two approaches are strong enough to lend credence to the existence of the transition. What seems 
less clear is whether the transition in the fluid involves the same conical structure as in [35^ [36] . 
The strong gravitational interaction between the images of localized black holes close to the merger, 
compared to the absence of interaction between the images of fluid balls, indicates the possibility 
of at least differences of detail in the two cases. 

An even more uncanny similarity in the phase structure of both fluid tubes and black strings 
refers to the critical dimension of spacetime for which the transition between uniform and non- 
uniform phases changes from first to second orderly] For fluid tubes ref. [20] found (and we have 
reproduced) that this happens at = 11, while for black strings it happens at D^, = 14 for the 
microcanonical ensemble [37] and = 13 for the canonical ensemble [38] ■ If we assume that, as in 
the AdS context, the d-dimensional fluids are to be related to D = d + 2-dimensional black holes, 
then the apparent closeness in the values of the critical dimensions is truly startling. 

The qualitative similarities between the dispersion relations in the Rayleigh-Plateau and in 
the Gregory-Laflamme instability have also been pointed out before for the case of static tubes 
and strings [1]. The parallels extend to further dynamical details, such as the fact that for static 

^In this and in the next section we shall use d and D, respectively, for the spacetime dimensions in which the fluid 
lumps and the black holes live. 
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solutions, non-axisymmetric modes (m 7^ 0) are stable in the plasma, while the GL instability also 
only affects s-modes [29l[39]. Furthermore, a boost along the static plasma tube has the same effect 
as boosting the black string: it is only a kinematic effect increasing the threshold wavenumber for 
s-modes and leaving all other modes stable [29j. 

Rotating black strings are affected by the same instability. The GL instability in rotating black 
strings was analysed in [JO]. They concluded that the GL instability persists for rotating black 
strings all the way to extremality and the threshold wavenumber increases as the rotation grows 
(see Fig. 1 of [40j). We could only treat small rotation analytically, but our dual plasma results, 
summarized in Fig. [U agree with this: the instability persists as rotation grows and it actually gets 
stronger, with the critical wavenumber increasing with the rotation. In fact this can be understood 
intuitively. The centrifugal force scales with the radial distance, so it is larger at the crest of the 
non-uniformity triggered by the RP instability than at the trough. As a result the distance between 
the crest and trough is stretched: rotation enhances the instability. An interesting qualitatively 
new feature that rotation adds is that non-axisymmetric modes (m 7^ 0) can become unstable to 
the Rayleigh-Plateau instability for sufficiently large rotation. 

Among the phases for rotating tubes that we have uncovered, a particularly interesting one is 
the pinched non-uniform tube (pNUT), since it has no analogue for static tubes. Its shape certainly 
suggests that we should think of it as the result of a bulge in the non-uniform tube beginning to 
behave like a rotating ball that develops a pinch. Since pinched black holes are also expected 
to exist in six-dimensional vacuum gravity (actually in any dimension larger than five) [321 [33], 
it is then natural to conjecture that rotating black strings in six dimensions also have phases of 
pinched non-uniform black strings, in addition to the already known uniform black strings, uniform 
black tubes {i.e., the black ring of [3T] times a flat direction), and the expected more conventional 
non-uniform black string. 

Such remarkable similitude between the dynamics of fluid tubes and of vacuum black strings 
makes it tempting to try to push it further in order to address an outstanding open problem: the 
time evolution of the GL instability and its final fate. The final stage of this evolution is not known 
in gravity, since the available numerical code, which has been developed only for D = 5, breaks down 
for late times [H]. On the other hand, the time evolution of the (non-relativistic) RP instability has 
been the object of several experimental and numerical studies, see e.g., ^42j which considers viscous 
fluids in four spacetime dimensions (unfortunately, higher-dimensional fluids remain unavailable 
experimentally). We would expect that for static initial data, the classical and relativistic evolution 
should not differ substantially. It is then justified to compare the initial time evolution of the GL 
instability. Fig. 4 of [H], with the early stage dynamical evolution of the RP instability, e.g.. Fig. 1 
of [42j. There appears to be a nice match between the two. In the fluid, where the evolution has 
been followed until its endpoint, the uniform tube pinches-off and breaks. Starting from a single 
sinusoidal perturbation in a cylindrical liquid bridge (fluid tube), the higher harmonics generated 
by non-linear effects are responsible for the development of a long neck that breaks the tube in a 
self-similar process. One ends up with an array of satellite drops with different sizes [420. 




Actually, two main regimes have been observed depending on the viscosity of the fluid. For low viscosity, there 
is repeated stretching and breakup, ending with several satellite droplets. For high viscosity, a single breakup that 
takes a long time can occur and in the end a single satellite drop is formed. 
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We may expect such behavior to be representative of the time evolution of the instabihty of 
fluids in subcritical dimensions d < d=K = 11. A numerical study of the non-relativistic evolution for 
inviscid fluids in several dimensions has been made in [21]. The endpoint of the RP instability does 
appear to differ for subcritical and supercritical dimensions. For d < the endpoint is a drop or 
an array of drops. Again, it is tempting to argue in favor of a similar evolution in the gravitational 
system, in which the black string pinches off completely and forms an array of black holes, possibly 
of different size The same caveats apply here as in the discussion of the ball-tube/black hole- 
black string transition: the fluid description breaks down near pinch-off; and blobs of fluid do not 
self-interact whereas gravitating horizon blobs do. However, it is already significant that pinch off 
does occur for viscous fluids, and that no study reveals any sign of the instability slowing down as 
pinch-off is approached in d < . 

The situation may be cleaner above the critical dimension. Ref. [;21j finds that in this case 
the endpoint of the RP instability is a non-uniform tube with constant mean curvature. In this 
case there is no concern about the validity of the fluid description, since the tube can remain thick 
enough during its evolution. The role of viscosity is not too clear, but it might be expected to 
smooth the evolution. This suggests that the endpoint of the GL instability on a black string 
above the critical dimension could be a non-uniform black string, like ref. [43] suggested. 

The non-linear GL evolution when rotation is present has not been studied yet but again, the 
(non-relativistic) RP evolution of fluids is well studied (see e.g., |44j ) : rotation introduces twisting 
effects in the pinch off of the tube and as rotation increases a dramatic centrifugal ejection of drops 
is observed (see the photographs of [44J). 

Of course it is possible that black strings and black holes in vacuum gravity behave in a closely 
similar way to fluids only as long as one considers configurations that are stationary, slowly evolving, 
and/or not too non-uniform. The dynamics of vacuum gravity may depart significantly from fluid 
dynamics away from these regimes and furthermore, in contrast to the situation for gravity in 
AdS, it is not obvious whether it can be accurately pictured as fluid dynamics in some controlled 
approximation. The next section addresses this issue in more detail. 

8 The black hole/fluid analogy revisited 

The similarities discussed in the previous section naturally prompt the question of to what extent we 
can trust fluid dynamics as a guide for unknown dynamics of vacuum black holes in asymptotically 
flat space. In other words, since the analogy between fluids and black holes was observed in the first 
place for asymptotically flat black holes, is there any precise meaning to it? Clearly, we cannot just 
take the results for AdS black holes in the limit in which the cosmological constant vanishes. The 
fluid description for AdS black holes requires that their size be much larger than the AdS curvature 
radius, so black holes smaller than this size fall outside the scope of the dual hydrodynamics. 

We have discussed a number of qualitative or semi-quantitative features that fluid lumps share 
with black holes in vacuum gravity. However, there is at least one crucial difference that would 

These have already been constructed in [22] , It is entropically favorable for these arrays of black holes to merge 
into a single one. 
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appear to preclude a precise correspondence. Classical General Relativity in vacuum is a scale- 
invariant theory, which implies that all properties of a black hole scale uniformly with mass. That 
is, all Schwarzschild black holes are essentially the same; all Kerr black holes with the same value 
of J/GM^ also have the same properties; and all black strings in Kaluza-Klein space are the same 
when the compact radius is scaled to keep L/(GM)^/(^~^) constant, e.g., the Gregory-Laflamme 
threshold mode wavelength for a static black string must scale like the inverse of the horizon 
transverse radius. This scale invariance is not present for the fluid: the surface tension sets a 
scale that distinguishes fluid lumps of different size: a fluid ball of a given radius is not simply 
a scaled-up version of a ball of half that radius — for instance, the former can more easily break 
up than the latter. For a large plasma ball, the relative entropy cost of breaking it in two pieces 
becomes arbitrarily small as the radius R of the ball gets arbitrarily large. For instance, for a ball 
of fluid with initial entropy Si, temperature Tj and radius Ri, that breaks into two equal balls, each 
with entropy Sf/2, temperature Tj and radius Rf, keeping the total energy unchanged, the ratio 
between final and initial entropies is 



Si 



+ (8-38) 



fluid 



Tf \ p \Ri Rf 



where we have assumed that the ball radii are much larger than the length scale a/p set by the 
energy density p and surface tension a of the fluic0- To leading order, RJ^ = R^^ > R^^, 
and moreover, if the balls have negative specific heat the smaller final balls will be hotter than the 
initial one, Tf > Tj, with Tf/Ti — 1 a positive and small quantity of order a/{pRi). Then Sf < S,, 
so the process is suppressed, but the relative entropy cost of breaking 
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Sf — Si 



fluid 



^ -dd-^ (8.39) 

fluid 



is small, with large fiuid balls breaking more easily than smaller ones. Here 5d is a positive, 
dimension-dependent number which, for a plasma with the equation of state considered in sec. 12.21 
is 

5rf = (d-l)(2^-l) . (8.40) 

For a Schwarzschild black hole, in contrast, the relative entropy cost for splitting into two equal 
black holes keeping the total mass fixed, remains constant independently of the black hole size, 

_ Sf-^ 
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1 



= 2"^-l. (8.41) 

bh 



It is amusing that, if we make the identification D = d + 2 that is suggested by the duality between 
SS-AdSd+2 black holes and d-dimensional fluid lumps, the powers of 2 in these formulas coincide. 
But on the other hand, the identification between (j8.39p and (j8.4ip appears to require that a ~ pR. 
This is not only a rather bizarre behavior for the surface tension, which grows linearly with the 
size of the fluid ball. It also implies that the black hole is equated with a fluid ball that is always 
smal|^, in the sense that surface effects are always strong so the regime (|6.34|) is never attained. 

^^To leading order we can neglect the difference in energy density in the initial and final fluid balls. 
""^^This is indeed the case for small black holes in AdS, which resemble asymptotically flat ones. 
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The apparent reason for this is, again, that vacuum gravity does not have a scale parameter to 
characterize a black hole as parametrically large. 

So the crux of the problem is that we are trying to relate objects in parameter spaces of different 

dimensionality! For the black holes, as argued, we are always free to fix a scale so as to set the 

1 

mass, or the length {GM) 0-3 ^ equal to one, and then completely characterize the black hole by 
its angular momenta. For a fluid lump, the theory comes with a scale already, namely the length 
scale a/ p, so in order to specify a fluid ball we must provide its size relative to that scale (or its 
energy), and its angular momenta — i.e., one more parameter than for a black hole. 

One may point out that a natural scale that appears in black hole physics is the Planck length, 
which does allow to make a distinction between small and large black holes. However, this is of no 
help to make a better comparison between vacuum black holes and fluid lumps, since in the limit in 
which the 'cutoff' scale becomes negligible, namely, when the fluid lumps or the black holes become 
very large relative to that scale, their properties are clearly different, as illustrated above. It would 
be desirable to have a way of deciding whether a black hole is large or small within classical vacuum 
General Relativity, but this seems hard in the absence of a classical fundamental scale. 

Does this, then, imply, that there can be no limit in which the equations of the dynamics and 
the phase space of vacuum black holes can be mapped into the equations and phase space of a fluid, 
with effective parameters that satisfy 

(8.42) 

In the following we make some observations that suggest otherwise. 

Let us first note that an obvious difference between fluids and black holes is that two discon- 
nected lumps of fluid do not attract each other (at least within the framework in which we are 
considering fluids) . This must obviously imply that some phenomena, like the deep non-linear evo- 
lution of the GL instability, or the mergers involving vacuum black holes or black rings, may behave 
differently than those for their fluid counterparts: in the former case the gravitational attraction 
between different lumps on the horizon presumably plays a role similar to the attraction between 
two black holes and thus affects their evolutioiJHI- Thus, any limit in which vacuum black holes 
behave accurately as fluids must be a limit in which this gravitational attraction is suppressed. 

Such a limit seems to arise if we consider gravity in a spacetime with a large number of non- 
compact dimensions. As the number of dimensions increases, the gravitational potential becomes 
steeper and more localized near the source, and flatter and weaker at larger distances. Furthermore, 
since the black hole entropy 

S oc (8.43) 

is directly proportional to the mass in the limit D — > 00, there is no entropy cost in splitting a 
black hole: (I8.4ip yields AS/S — > in this limit. As we observed above, this is just like for a large 
drop of fluid. Indeed, both for a fluid with vanishingly small surface tension and for a black hole 
in D — > 00 the (Hagedorn-like) relation E = TS is satisfiecj^. 



^*That these effects must, of necessity, be strongly suppressed for the class of AdS black holes that are dual to 
plasma lumps, is a remarkable consequence of the duality. 

^^For the fluid lump E here is the difference between its energy and the energy of the (unconflned) vacuum it 
displaces. 
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This enhanced instability of the black hole horizon to break up has a rather precise manifestation 
in the equivalence between the Gregory-Laflamme and Rayleigh-Plateau instabilities. Given the 
observations above, we can expect a black string in D ^ oo to be able to split into fragments of 
any length, i.e., the wavelength of the threshold mode is expected to approach zero. Indeed, the 
GL instability in arbitrary dimension has been studied in [45] . and even if the equation for the 
threshold mode is fairly complicated, it simplifies greatly if we expand in 1/D to leading order, 

X"ir) + -x'{r)-k\{r)=0. (8.44) 
r 

This is essentially the same as the equation for the RP threshold (w = 0) mode (j5.2ip . and it 
becomes precisely the same if we identify ~ n ~ c? for large D and d. The wavelength of the GL 
threshold mode in this limit vanishes as r^/V^, in accord with the expectation above. 

Expanding on these observations, we can try to identify what plays the role of the dimensionless 
quantity a/{pR) in the gravitational side. The answer is obvious: since the only small dimensionless 
parameter we have at our disposal in the gravity side is 1/D, both must be related. A more precise 
mapping can be obtained by comparing (j8.39p and (|8.4ip . For large D, the latter becomes 



AS 
~~S' 



(8.45) 



which leads to the identification 



We are not necessarily proposing that the limit d ^ oo must be taken in the fluid side too, but if 
we do, and use the value (|8.40p . 6^ log 2, then the identification becomes precise, ^ = D^^, at 
least to leading order in 1/D. 

Eq. (I8.46P illustrates a main feature of our proposal. As we have discussed, the main problem 
for a fluid/black hole correspondence is the apparent absence in classical vacuum GR of a way to 
decide whether a black hole is intrinsically large or small. Our suggestion amounts to saying that 
a black hole is large or small depending on the number of spacetime dimensions it lives in. Black 
holes in low dimensions are to be regarded as tiny droplets, with surface energy comparable to 
bulk energy, whereas black holes in very high dimensions behave rather like large lumps of fluid. 
In other words, we are using D as an additional parameter characterizing a black hole. We are 
making crucial use here of the observation that in classical General Relativity in vacuum the only 
dimensionless tunable parameter at our disposal is D, a point that has also been urged by B. Kol 
and collaborators [36] . 

It might be that these similarities between fluids and black holes at large D only reflect purely 
geometrical aspects of the relation between volumes and areas at very high dimensions (note in 
particular that D, but not a/{pR), appears in (j8.44p ). But at any rate we feel that these observa- 
tions suggest that perhaps by considering a 1/D expansion, and possibly accounting for the effects 
of gravitational interaction in this expansion, a more accurate map betwen vacuum black holes and 
fluids may be achieved that allows to understand not only the similarities that have been observed, 
but also the differences. 
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